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Abstract To achieve a maximum compression of system matrix in positron emission tomography (PET) image 
reconstruction, we proposed a polygonal image pixel division strategy in accordance with rotationally symmetric PET 
geometry. Geometrical definition and indexing rule for polygonal pixels were established. Image conversion from 
polygonal pixel structure to conventional rectangular pixel structure was implemented using a conversion matrix. A set 
of test images were analytically defined in polygonal pixel structure, converted to conventional rectangular pixel 
based images, and correctly displayed which verified the correctness of the image definition, conversion description 
and conversion of polygonal pixel structure. A compressed system matrix for PET image recon was generated by tap 
model and tested by forward-projecting three different distributions of radioactive sources to the sinogram domain and 
comparing them with theoretical predictions. On a practical small animal PET scanner, a compress ratio of 12.6:1 of 
the system matrix size was achieved with the polygonal pixel structure, comparing with the conventional rectangular 
pixel based tap-mode one. OS-EM iterative image reconstruction algorithms with the polygonal and conventional 
Cartesian pixel grid were developed. A hot rod phantom was detected and reconstructed based on these two grids with 
reasonable time cost. Image resolution of reconstructed images was both 1.35 mm. We conclude that it is feasible to 


reconstruct and display images in a polygonal image pixel structure based on a compressed system matrix in PET 


image reconstruction. 
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1 Introduction 


Nowadays, as an essential branch of molecular 
imaging technologies, positron emission tomography 
(PET) has made rapid developments. Based on modern 
high-sensitivity 3D data acquisition mode and 
statistical iterative reconstruction algorithms, 
high-quality molecular and functional images from 
organs to cells of human bodies can be obtained. PET 
plays an important role in oncology, cardiology and 
neurology diagnosis, including early detection, stage 
evaluation, and 


of malignant post-treatment 


monitoring of human diseases. 


It remains a great challenge though to 
fully 3D PET 
reconstruction algorithm in a clinical PET system, 


implement a high-resolution, 
since there are numerous number (about 10^)! of 
lines-of-response (LORs) events need to be calculated 
when a modern PET machine works in 3D mode. 
Correspondingly, a huge-sized (about 10? elements) 
system matrix is needed in image reconstruction. In 
order to address this issue, it is critical to compress the 
system matrix and obtain best trade-off between image 
quality and reconstruction speed. 

Considering that most modem clinical PET 
systems employ a symmetric ring-shape geometry, a 
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concept of taking advantage of this symmetry with 
irregular image pixel shape instead of conventional 
rectangular pixels to achieve a maximum compression 
of system matrix was put forward, which recently 
becomes a new focus area in fully 3D PET iterative 
algorithm development. Compared with other methods 


è ; 2 
of system matrix compression"! 


, compression with 
symmetric is direct and has been proven to be efficient 
to achieve better compression factor. 

Early approaches designed a polar grid based 
system matrix structure which has rotational and 
reflectional symmetries. For example, such principle 
was proposed to accelerate PET iteration process by 
polar pixel structure), 

Recent years, with the popularization of 3D 
PET technology and explosive growth of acquisition 
data size, more attention was involved in this field. In 
addition to the polar grid, various irregular pixel 
shapes were presented in commercial or laboratory 
PET systems such as in MicroPET Focus 120°! and in 
LabPET!! A software toolkit PRESTO!!! was 
developed, which integrated the generic calculating, 
storing and accessing techniques based on the 
rotational symmetry in a C++ library to improve the 
efficiency and accuracy of image reconstruction. 
Comparing to the original polar pixel based geometric 
structure, an irregular division can overcome the 
unevenness of pixel area distribution across the image 
FOV. 

It has been demonstrated that a polar voxel 
with square-like shape produces fewer artifacts ^ 9], 
so one usually has to design specific division mode for 
each compression factor as shown in Ref.[11]. The 
increased complicity of pixel shapes in the irregular 
division mode brings further challenges in image 
reconstruction and display. 

Owing to that most existing PET systems have 
a polygon ring shape rather than a circle, a polygonal 
pixel grid may have an edge than traditional polar grid. 

This work proposes a polygonal image pixel 
division strategy based on the characteristics of 
rotationally symmetric PET geometry. It is adaptable 
to the structure symmetries of various types of PET 
systems. One can simplify the calculation of the 
system matrix at a compressed format by the usage of 
rotational symmetry of pixel structure during the PET 


image formation process. The proposed method was 
applied to a practical small animal-PET system. 

In view of the definition of this pixel grid and 
calculation of the system matrix, we accomplished the 
iterative reconstruction and image display of 
experiment data of a hot rod phantom with 30 mm 
both in diameter and length based on the practical PET 
scanner. 3D detected data were organized into 
sinograms. System matrix was calculated by tap model 
and compressed by rotational and axial symmetries. A 
Cartesian grid based image reconstruction algorithms 
with same pixel size was also generated with 
tap-model. Storage size of system matrix, time cost 
and image quality were compared for these two 


methods. 


2 Methods 


2.1 Rotationally symmetric polygonal pixel based 
image expression & display 

The core concept of this paper is based on a novel 

grid. A 

representative division pattern of the polygonal pixel 


rotationally symmetric polygonal pixel 
grid is shown in Fig.l(a) The structure contains 8 
segments in Fig.1(a) which are rotationally symmetric 
with respect to the center of the image. The 
arrangement of polygonal pixel is same from segment 
to segment. 

Compared to the conventional rectangle pixel 
division pattern (Fig.1(b)), the 
rotational symmetry helps to compress the storage 


introduction of 


space of the system matrix. 

Compared to the existing polar pixel grid 
pattern (Fig.1(c)), the regularity of this division mode 
ensures a simplified expression and automatic 
segmentation. In another aspect, there are only three 
types of basic pixel shapes. One can easily adjust the 
area of each type of pixel to minimize the difference of 
pixel area over the image FOV. Therefore, this 
structure may have an edge over polar grid in terms of 
reducing pixel area unevenness. 

Enslaved by the typical image display mode of 
modern computer, the main challenge of this novel 
structure is that the image expressed in this rotational 
symmetry polygon pixel grid has to be translated into 


the cubic ones for image display. Besides, a new 
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method is needed to calculate geometrical weights of 
system matrix since the universal algorithms such as 
Siddons ray-tracing method!" is not applicable in this 


particular case. 


2.2 Geometrical description of polygonal pixels 


Based on the concept of building arotationally 


symmetric geometrical structure, the definition of 
polygonal pixels consists of three steps: (1) segment 
definition, (2) layer definition inside segments and (3) 
pixel definition, as illustrated in Figs.2(a), (b) and (c), 
respectively. 


Fig.1 Comparison of basic structures of (a)rotationally symmetric polygonal pixel, (b)conventional rectangle pixel and (c) polar 
pixel. The compression factors of polygonal and polar pixels are both eight. 
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Fig.2 Division mode of rotationally symmetric polygonal pixel grid. (a) segment-division. P -rays starting from the image center 
(origin) define the P fan-shaped segments (0~P-1, in (a), P=8). Angle of two adjacent rays is also the vertex of each segment defined 
as 0. An index number p = 0, 1... 7 is anti-clockwise assigned to each of the 8 segments. The starting and ending boundaries of only 
segment 0 and 7are shown for clarity. (b) layer-division. A group of parallel lines which are perpendicular to the vertex angle bisector 
divide a segment into several layers. The layer-height p is defined as the length of intersecting lines of adjacent layer-boundaries and 
the segment-boundary. The indexing number of layers is increased from the layer containing the vertex. (c) pixel-division. A group of 
lay-intersecting lines which are parallel to the vertex angle bisector divide pixels in-layers. The layer 3 has 4 pixels. Inside layer 3, the 
0-th pixel contains the ending boundary, and the 3-th pixel contains the starting boundary. 


Segment definition: When the compression 
factor is P, P rays are drawn starting from the center of 
image field of view (FOV), which divides the FOV 
into fan-shaped segments. The angle 0 between the 
two adjacent rays, ie. the vertex angle of each 
segment is @=2n/P. The shape and area of every 
arrangement is equal to each other. 

Starting from the positive x-axis which is 
defined to be the starting boundary of the first segment, 
an index number p-0, 1...P-l is anticlockwise 


assigned to each of the P segments. The ending 
boundary of p-th segment is the same as the starting 
boundary of (p-1)-th segment, where p=0, 1...P-2. 
The ending boundary of segment P-1 coincides with 
the positive x-axis. 

As all the segments are identical in shape, the 
arrangement of pixels can be defined only once for 
one segment, which can be rotationally replicated to 
other segments. Therefore, the definition of layers and 
pixels are discussed only in segment 0. 
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Fig.3 In the segment belonging to a 8-segment polygonal 
pixel grid, of which the vertex angle is 0, the layer-height is p 
and the layer number is S, there are only three shapes of 
polygonal pixels. There are only one pixel shaped like triangle 
and the area is j!sinO/2 . There are 2S right angle 
trapezoidals, of which the area is. 3p? sin / 4, and the number 
of rectangle is (s-1)s/2 as well as an area of p sinÓ. 


Layer definition: The area inside one segment is 
divided into 
layer-boundary lines which are perpendicular to the 


several regions by a group of 


vertex angle bisector and intersect the 
segment-boundary with equal distance p. The distance 
between two adjacent  layer-boundary lines 
is p cos(0 / 2) . The area embraced by the two adjacent 
parallel lines and the two segment-boundary rays is 
called a layer, except the centermost layer which is 
formed by the two segment-boundary rays and one 
layer-boundary line which is the closest to the vertex 
angle. Each layer is identified by an index number s. 
The layer containing the vertex is layer 0, whose 
adjacent layer is layer 1, and so on. The total number 
of layers is defined as S, which is determined by the 


size of the image. 


Table 1 Parameters of three basic geometric structures belonging to an S-layer segment of the rotationally symmetric pixel grid 


LN 


m 


| | 


Parameters 

Shape name Triangle Right angle trapezoidal Rectangle 
Number of pixels 1 28 (S-1)S/2 
Area p sin@/2 3p! sin@/ 4 p sinb 
Relative area ratio 2 3 4 


Pixel definition: As illustrated in Fig.2(c), a set 
of lines parallel to the vertex angle bisector are drawn 
which divide each layer into one or multiple polygonal 
regions. Those lines are denoted as lay-intersecting 
lines to be distinguished from layer-boundary line in 
layer definition, and each polygonal region is defined 
as one pixel. The distances between the neighborhood 
layer-intersecting lines are all equal to 2psin0/2. 
The s-th layer has s+1 pixels that are in a row. The 
row-index for those pixels labeled as ¢. We define that 
in each layer, the 0-th pixel contains the ending 
boundary, and the s-th pixel contains the starting 
boundary. 

With the above steps described, one can 
establish a rotationally symmetric polygonal pixel grid 
and enumerate all the pixels with three parameters, i.e. 
segment-index p, layer-index s and row-index t. 
Therefore, a pixel is uniquely specified by (p,s. t). 

In comparison with conventional 
rectangular-pixel based image division pattern, the 


proposed pixel division strategy has a P-fold rotational 
symmetry when there are P segments, leading to a 1/P 
storage size of the system matrix. The existing 
polar-pixel based methods can achieve a similar 
compression factor. However, considerable difference 
of both pixel shape and area exists in pixel division 
pattern. E.g. in Ref.[15], the largest area ratio of pixels 
reaches almost 5, which may lead to non-uniform 
image quality distribution and even image artifacts. As 
listed in Table 1, in the proposed polygon pixel 
division strategy, there are only 3 basic pixel types: 
triangle, right angle trapezoidal and rectangle, and the 
area ratios of these three types of pixel are 2:3:4. The 
simple pixel shapes and reduced pixel area difference 
are beneficial in minimizing resolution non-uniformity 
and reducing image artifacts. 


23 Pixel indexing and addressing 


In order to store a digital image into computer memory, 
a pixel indexing rule is needed to map each pixel to a 
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memory unit. Likewise, a pixel addressing method is 
needed to calculate the pixel coordinate given its 
memory address. 

In an image expressed by rotationally 
symmetric polygonal pixels, each pixel is uniquely 
represented by its three indices (p,s,f). When the image 
is stored in computer, it occupies a continuous space of 
computer memory, for each pixel (p,s,f) a unique 
address i needs to be determined. 

Figure 4 shows a representative scheme for 
assigning memory address to image pixels in an 
8-segment, 5-layer image. The pixel indices are 
determined with a  segment-by-segment and 
layer-by-layer sequence. The pixel indices inside p-th 
segment are from 15p to 15p+14. Inside a segment, the 
pixel indices are increasingly assigned from smaller 
layers to larger layers and from smaller rows to larger 


rows inside each layer. 


Fig.4 Diagram for stored orders (pixel code) of pixel values 
of the rotationally symmetric polygonal structure, taking P=8, 
S=5, F-120 as an example. 

According to above rules, one has the mapping 
relationship between i and (p,s, t) is: 

s(s +1) 
2 
Where J=(S+1)(S+2)/2 is the total number of 

pixels in one segment. 
Meanwhile, the coordinate (p,s,f) of the i-th 


pixel is determined as: 
i 
= — 2 
p l ;| (2) 


i=Ix pt 


+t (1) 


"e Bees sae?) a 


_ s(s +1) 


t-i-pxl (4) 


2.4. Image pixel structure conversion 


To view an image stored in the rotationally symmetric 
polygonal pixel structure on computer screen, it has to 
be firstly converted to a conventional rectangular pixel 
structure. The conversion process can be expressed as: 


Image,, = Q,,,Image, (5) 


Where Image,, is the converted image in rectangular 
pixel based structure, Image, is the original image in 
rotationally symmetric polygonal pixel structure, both 
organized in column vector format, and Qyy is the 


conversion matrix. 


Fig.5 Sketch map for pixel transforming. The m-th square 
pixel, whose area is Sm, is covered by several pixels in 
rotationally symmetric structure; the value of which can be 
expressed as a weighted sum of those unconventional pixels 
covering this square one. The weight is considered as the ratio 
of the overlap area and the total area of the cubic pixel: 
Qmi^ Sinn! Sm 

As illustrated in Fig.5, the value of m-th 
rectangular pixel in Image,, can be expressed as a 
weighted sum of the values of multiple polygonal 
pixels overlapping it. The overlap area of the n-th 
pixel in polygonal structure and the m-th pixel in 


rectangular structure is denoted as $,,. The weight is 
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considered as the ratio of the intersection area and the 
total area of the cubic pixel. Therefore, one has 


N-1 

Ys. ^ d, 
b = n=0 — mn 
m S, S, 


2.5 Verification of image expression and 


structure conversion 


In order to verify that an image is correctly stored in a 
rotationally symmetric polygonal structure and 
converted to a rectangular pixel structure, a set of test 
images are analytically defined: 


a=0; s-0, 2, 4 


(a) 


a(p,s,t)=0, where s is odd 


=], where s is even (7) 
By selecting different segment number, a set of 
images A, is numerically defined. Fig.6 illustrates the 
test images of 8, 9 and 12 segments, respectively. 
Those images are converted to B,, stored in 
rectangular pixel structure, and displayed on computer 
screen for examination. If the structure is defined 
correctly and the transforming matrix works 
effectively, one should see the images corresponding 
to the predefined ones in various rotationally 
symmetric polygonal grids. 


(b) (c) 


Fig.6 Schematic diagram for the testing images of 8 segments(a), 9 segments(b) and 12 segments(c) for image expression and 
conversion, where the value of pixels equals 0 when s=0,2,4; and the value of pixels equals 1 when s=1, 3. 


3 Generation and verification of system 
matrix 


3.1 Generation of a compressed system matrix 


In the calculation of the system matrix, each 
line-of-response (LOR) is considered as a tape with 
certain length. As shown in Fig.7, each LOR is 
uniquely identified with its projection angle g and the 
distance to the center of FOV r. Therefore, the j-th 
LOR is denoted as j(g,r). Each element c;; in the 
system matrix records the probability that the j-th LOR 
detects a photon emitted from the i-th pixel, denoted as 
i(p,s,t). In this study, the impact of physical effects, 
such as positron range, photon non-collinearity, 
attenuation and scattering are ignored, and pure 
geometric detector response is modeled. In this way, 
cijis considered to be proportional to the overlapping 
area of the i-th pixel and j-th LOR. 


It is clear the rotational symmetric structure in 
the image domain can be utilized to compress the 
system matrix. As shown in Fig.8, considering two 
image pixels ig and 7, when they have same layer index 
s and row index f, and when the angle from LOR joto j 
is equal to pO, i.e. the angle from segments poto p, the 
overlapped area between pixel i and LOR j is 
rotationally symmetric to that between pixel ip and 
LOR jo, thus the corresponding system matrix element 
has same values, ie. c =C} 

The relationship between i and ig is shown 
followed: 

i, =i-1xp (8) 

The relationship between j and jọ can be 
expressed by relationships of r and rọ as well 


as ọ and g,: 
har (9) 
—J9-pO  o-p070 
P = { Q-pO*z g-p0<0 (10) 
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Fig.7 General view for calculation of system matrix with projection-tapes model. In PET detector structure, The ratio of the 
projection tape ji.) and i-th pixel’s overlapping area Sli, and the total area of pixel i, S;, are taken as the weight of system matrix 
when we regarded projection lines as a group of tapes with a certain width. 


3.2 Verification of system matrix 


The following verification studies are based on a 
simulated typical small animal PET system. As shown 
in Fig.9, the detector consists of 6 rings, on each of 
which 32 detector blocks are evenly distributed. Every 
block is composed of 8x8 LYSO crystals. The size of 
crystal is 1.59 mmx1.59 mmx10 mm. The detector 
ring diameter is 147.2 mm, and the size of FOV is 
30 mm in diameter and 76.32 mm in axial length. 
Radial center of the entire detector field of view is a 
circular with a radius of about 30 mm. 

A rotationally symmetric polygonal pixel 
structure with P=12, p=0.5 mm is defined. The 
layer-height p is chosen so that it is much smaller than 
the desired PET system resolution i.e. 1.6 mm. To 
cover the 30 mm FOV, the image has 60 layers. The 
acquired 2D PET data are organized to a sinogram 
with 155 angular samplings over 180 degree and 60 
radial samplings over 30 mm trans-axial FOV, and the 
width of the LOR tape is the same as the layer-height. 


By defining different source distributions in the 
image structure, the corresponding sinogram data were 
generated by projecting the source image into LOR 
space using the pre-calculated system matrix. The 
sinograms are visually examined and compared to 
theoretical expectation to verify the correctness of the 


system matrix and projection model. 


4 Image reconstruction 


Data for image reconstruction were acquired with an 
experimental small animal PET system. The detector 
consists of 5 rings, on each of which 24 detector 
blocks 
composed of 9x9 LYSO crystals. The size of crystal is 


are evenly distributed. Every block is 


2.11 mmx2.11 mmx10 mm. The inner detector ring 
diameter is 76.33 mm, and the size of FOV is 50 mm 
in diameter and 95 mm in axial length. 

A rotationally symmetric polygonal pixel 
structure with P=12, p=0.33 mm was defined. The 
layer-height p was chosen so that it was much smaller 
than the desired PET system resolution. To cover the 
50 mm FOV, the image had 79 layers. 
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Fig.8 Predominant feature in symmetries can be utilized to calculate the weights of pixels in no-zero segments to all projections 


tapes. The weight c, of the i-th pixel (i>S) for projection tape j,,,, is same as the weight c 


when two image pixels ip and i 


iosjo ? 


have same layer index s and row index t, and the angle from LOR jo to j is equal to pO, i.e. the angle from segments py to p. 


ro 


Fig.9 Diagram for the detectors axial sectional view and the overall structure in a typical small-animals PET. 


The acquired 3D PET data were organized to a 
sinogram with 108 angular samplings over 180 degree 
and 70 radial samplings over 46.76 mm trans-axial 
FOV, and the width of the LOR tape was the same as 
the layer-height. 

After the establishment of the system matrix, 
the projection data acquired from the real PET scanner 
were reconstructed based on the rotationally 
symmetrical structure. Using the image transformation 
matrix defined in Eq.(5), the reconstructed images 
were transformed and displayed for visual 
examination. 

We also conducted the image reconstruction 


based on a Cartesian pixel grid with 142x142 pixels of 


0.33 mmx0.33 mm area. The system matrix was 
calculated using the same tap model for LORs. 
The commonly used OS-EM algorithm"! 
was chosen for PET image reconstruction in this study. 
The OS-EM algorithm is based on maximum 
likelihood estimation of Poisson distributed signal. 
The objective of the iterative algorithm is to maximize 
cost function which defines a logarithmic likelihood 
function between the projection of the image estimate 
and the actual projection values. An iterative 
expectation maximization algorithm is used to 
optimize the cost function. The projection views are 
grouped in L subsets, the image is updated after each 
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subset is considered. The iterative formula is as 
followed: 


Where p; is a simulated detected photon number in j-th 
LOR, f; is image value of the i-th pixel after the k-th 
iteration when using projection belonging to s-th 


subset, and c;; is an element in the system matrix 
records the probability that the j-th LOR detects a 
photon emitted from the i-th pixel. 


Fig.10 Display results of the test images defined in Eq.(7) converted to Cartesian coordinate. The segment number P of figure (a), 
(b), (c) is respectively 8, 9 and 12, and the pro-defined images are shown in Fig.6. The layer-height p of the rotationally symmetric 
structure is 3mm, and. The pixel size of the rectangular pixel structure is 1mm, and the final size of the image is 128x128. 


The experiments in this study were based on a 
practical multi-ring small animal-PET scanner. A 
Derenzo phantom of F-18 source with 6 groups of hot 
rods was scanned during experiments. The diameters 
of the hot rods were 0.75 mm, 1.0 mm, 1.35 mm, 1.7 
mm, 2.0 mm and 2.4 mm respectively. In each group, 
the distance between adjacent rods was twice of the 
rod diameter. List-mode coincidence event data were 
collected and histogrammed to a 2D sinogram 
projection. According to these data we can implement 
image reconstruction with calculated system matrix for 
both the polygonal pixel grid and cubic one in this 
study. Storage size of system matrix, time cost and 
image quality of can be compared for these two 
methods. 


5 Results and discussion 


5.1 Structural testing for rotationally symmetric 


polygonal pixel grid 


As defined in Eq.(7), three test images are defined 
with three different segment numbers: P=8, P=9 and 
P=12. These images are defined with rotationally 
symmetric polygonal structures, and converted to 
images expressed with conventional rectangular pixel 
structure as described in Eq.(5). The pixel size in the 


converted images is almost 1/5 of that of original 
images. Three converted images are shown in Fig.11. 
No image quality loss is observed after the conversion, 
which demonstrates the correctness of the image 
structure conversion process. 

Figure 10 displays results of the test images 
defined in Eq.(7) converted to Cartesian coordinate. The 
segment number P of Figs.10(a-c) 1s respectively 8, 9 and 
12, and the pro-defined images are shown in Fig.6. The 
layer-height p of the rotationally symmetric structure is 
3 mm. The pixel size of the rectangular pixel structure is 
1mm, and the final size of the image is 128x128. 


5.2 Validation of the generated system matrix 


through forward-projection process 


In this test, three different distributions of radioactive 

sources were defined: 

1) A point source located at the center of FOV; 

2) A point source located radially 29.6 mm off the 
center of FOV; 

3) A volumetric source which uniformly fills the 
entire FOV. 

Three corresponding sinograms were generated by 

forward projecting the above source images with the 

generated system matrix, which are shown in Fig.11. 


030202-9 


YU Yunhan et al. / Nuclear Science and Techniques 24 (2013) 030202 


180° 


g 90° 


09 


0 30 -30 0 30 


r/mm 


Fig.11 The results of the forward-projection test for System Matrix, (a) the sinogram for a point source located at the center of FOV, 
(b) the sinogram for a point source 29.6 mm off the center of FOV and (c) the sinogram for a volumetric source uniformly fills the 


entire FOV. 


In Fig.11(a), the sinogram is expected to be a 
straight line (r=0), and the generated sinogram agrees 
well with the expectation. When the point source is 
29.6 mm off-centered, the expected sinogram is a 
sinusoidal curve with 29.6 mm amplitude. In the 


direction of the generated sinogram is 29.5 mm. Since 
the size of FOV is 30 mm, the covered range well 
matches the FOV size. Therefore, the sinograms 
generated by system matrix are corresponding to the 
theoretical expectations, and it approves that the 


generated sinogram, the measured amplitude is system matrix is calculated correctly. 


| ry 
(b) 


Fig.12 Comparison of central transaxial and coronal slice of reconstructed phantom images based on a polygonal pixel structure 
with P=12 and p=0.5 mm (a) and a 141x141 cubic pixel grid (b). We implemented OS-EM iterative image reconstruction with 54 
iterations. Data were collected from actual list-mode data and organized into projection in 108 projection angles and 70 intervals. The 
sizes of derenzo model are 0.75 mm, 1.0 mm, 1.35 mm, 1.7 mm, 2.0 mm and 2.4 mm, respectively. 


29.6 mm. In Fig.11(c), the covered range in the r 


5.3 Image reconstruction results polygonal pixel grid and cubic grid are shown in 
Figs.12(a) and (b), respectively. 
A hot rod phantom with 6 sections of hot rods (0.75 : ; 
Figure 12(a) shows reconstructed images based 
mm, 1.0 mm, 1.35 mm, 1.7 mm, 2.0 mm and 2.4 mm . , , 
mM . on a rotationally symmetric polygonal structure, with 
in diameter, respectively) and 30 mm cross sectional ; 
. . . Á 12 segments and 79 layers. The layer height was p 
size was scanned using a practical animal PET system. : : 
im —0.33 mm. After the reconstruction, the images were 
2D sinograms (108 angular bins and 70 radial bins) : 1 
. converted to conventional rectangular pixel structure, 
were generated from the list-mode event data. The . : ' : 
: . . : . with 0.1 mmx0.1 mm pixel size. We can see typical 
images reconstructed with OS-EM iterations in 
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rotational polygonal pixel shaped artifacts in the lower, 
middle part of the zoomed image section. 

The right column shows reconstructed images 
based on a cubic pixel grid with 142x142 pixels of 
0.33 mmx0.33 mm sizes. 

Both in Figs.12(a) and (b), the 1.35 mm hot rod 
sections can be distinguished. Line profiles in Fig.13 
illustrated that the desired 1.35 mm image resolution 
can be achieved with the proposed rotationally 
symmetric polygonal structure. 


——— Polynoal Pixel Grid 
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Fig.13 Vertical line profiles through the section with 1.35 mm 
in diameter in the reconstructed phantom images in Fig. 12. 

Based on this small animal PET system 
mentioned previously, the memory size of system 
matrix can be compressed to 46.1 MB when elements 
are stored as 32-bit float numbers using this polygonal 
pixel grid. In comparison, when we use the 
conventional rectangular pixel grid, the required 
memory size reaches 581.5 MB. The compression 
ratio of the system matrix size is 12.6:1. 

It took 2804 s to run 54 iterations in polygonal 
pixel grid running in a PC with a 2.5 GHz Inter(R) 
Core(TM) i5-3210M CPU and 4 G memory. The 
iterations conducted in cubic grid cost 2440 s. 

Considering the total pixels number in 
polygonal gird is 1.88 times more than the regular one, 
the slight increase reconstruction time is acceptable. 


6 Conclusion 


In this work, we designed a novel rotationally 
symmetric polygonal pixel structure for PET image 


reconstruction. 


No image quality loss was observed after the 
conversion of image in polygonal pixel structure to a 
conventional rectangular pixel for image display when 
the area of rectangular pixels is 1/5 of the polygonal 
ones. Three different distributions of radioactive 
sources were forward-projected to sonograms and 
verified the correction of system matrix. 

Images were successfully reconstructed based 
on 3D list mode data sinograms of hot rod phantom 
from a practical small-animal PET. The image 
resolution of reconstructed phantom image is 1.35 
mm. 

A compress ratio of 12.6:1 of the system 
matrix size was achieved on this small animal PET, 
comparing with the conventional rectangular pixel 
based one. 

With acceptable time cost, the image resolution 
based on polygonal pixel structure is as good as the 
cubic pixel structure based one. 

In further studies, we will expand our study to 
the optimization of division mode of rotationally 
symmetric structure and seek solution for the 
unevenness of pixel shapes and areas. We can also 
the 3D 
calculation to establish foundations for practical 
applications in the full 3D PET. 


apply this structure to system matrix 
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